Statistically Locked-in Transport Through Periodic Potential Landscapes 
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Classical particles driven through periodically modulated potential energy landscapes are pre- 
dicted to follow a Devil's staircase hierarchy of commensurate trajectories depending on the orien- 
tation of the driving force. Recent experiments on colloidal spheres flowing through arrays of optical 
traps do indeed reveal such a hierarchy, but not with the predicted structure. The microscopic tra- 
jectories, moreover, appear to be random, with commensurability emerging only in a statistical 
sense. We introduce an idealized model for periodically modulated transport in the presence of 
randomness that captures both the structure and statistics of such statistically locked-in states. 



Objects driven through periodic potential energy sur- 
faces face a myriad of choices: either they follow the 
driving force or they become entrained along any of the 
commensurate directions through the landscape. Vari- 
ants of this problem appear in areas as diverse as driven 
charge density waves [l| , electronic energy states in two- 
dimensional electron gases , atom migration on crystal 
surfaces [3, chemical kinetics, and flux flow in type-II 
superconductors Q . Quite recently, this problem was in- 
vestigated using a monolayer of colloidal spheres in 
flowing water as a model system and a square array of 
holographic optical tweezers to provide the periodic 
potential energy surface. Depending on the array's ori- 
entation with respect to the driving force, the spheres 
were observed to trace out a Devil's staircase hierarchy 
of commensurate directions through the array, with par- 
ticular directions being preferentially selected over cer- 
tain ranges of orientations Trajectories deflected by 
a preference for commensurability are said to be kineti- 
cally locked in to the lattice. 

This anticipated result was accompanied by two 
surprises. In the first place, not all kinetically locked- 
in states were centered on simple commensurate direc- 
tions. Still more surprisingly, particles' microscopic tra- 
jectories in high-order locked-in states did not consist of 
sequences of commensurate jumps, but rather consisted 
of seemingly random lower-order hops whose combina- 
tion, nonetheless, were commensurate. The appearance 
of statistical rather than deterministic commensurabil- 
ity suggests an unexpected role for thermal randomness 
in structuring transport through periodic potentials, and 
has been dubbed statistical lock- in 0. 

This Letter presents a highly simplified model for 
statistically locked-in transport through mesoscopic po- 
tential energy landscapes that nonetheless accounts for 
the emergence of combination jumps and their statisti- 
cal commensurability. In particular this model reveals 
how the potential energy landscape's structure and ex- 
tent establish the discrete spectrum of travel directions 
selected by biased random walkers. While our discus- 
sion is directed toward the purely classical behavior of 
flowing colloids, similar results should emerge for biased 
quantum mechanical hopping through arrays of potential 



FIG. 1: Schematic diagram of kinetically locked-in transport 
through a finite domain. Potential wells of radius a are ar- 
ranged on a square grid of lattice constant b. An applied force 
drives objects through the array at angle 9, but, because of 
their interactions with the wells, the particles are most likely 
to follow a commensurate direction '(/'(S) through the array. 
Deflection by thermal forces or quenched random scattering 
centers causes the trajectories to deviate from commensura- 
bility, which is recovered only in a statistical sense. 



wells or barriers. Finally, we argue that the kinetically 
locked-in state selected at a given orientation can depend 
sensitively on particle size. Statistical lock-in therefore 
presents new opportunities for continuously sorting het- 
erogeneous materials into multiple fractions simultane- 
ously. 

We first consider an array of optical traps, as shown 
in Fig. n The centers of the traps are on a square lat- 
tice with lattice constant b. The array is taken to extend 
indefinitely in the y direction, while its extent in the x 
direction is N lattice constants. Rather than attempting 
to account for the detailed structure of a practical trap 
array, we instead focus on the simplest model that cap- 
tures the experimental phenomenology. For this reason, 
we model the traps as circular regions of radius a < b/2 
centered on grid points. Once a particle enters a trap's 
domain, it is translated directly to the center. A particle 
passing outside the circle is assumed to be unaffected. A 
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"trapped" particle eventually escapes and is carried along 
to its next encounter by the driving force. This is con- 
sistent with the observed dynamics of colloidal particles 
flowing through holographic optical trap arrays j5(. 

In such a situation there are two distinct regimes we 
may consider. The first is the "ballistic" regime where 
most particles pass through the lattice without encoun- 
tering traps. This corresponds to the limit of small trap 
density, rj = a/b. For sufficiently large values of 77, traps 
are packed closely enough together that particles' trajec- 
tories consist mostly of jumps from one trap to another. 
This we term the "lattice gas" regime. We shall first de- 
rive the condition on 77 for being in the lattice gas regime 
and consider only the case that meets this condition for 
the rest of the discussion. Consider the array as shown 
in Fig. n Consider a line drawn from one of the lattice 
points on the left edge which we will take to be the origin. 
Let the line have a slope m. Now the shortest distance, 
s, between the line and a generic point on the lattice, 
represented by (p, g), with p G {1, • • • , N} and g G Z, is 
given by 



b\q — pm\ 
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(1) 



We now wish to compute the largest possible value, Sm, 
of the distance from the line to the lattice point nearest 
to the line for all possible slopes. If the trap radius a 
is larger than this distance then all straight line trajec- 
tories originating at a lattice point on the left side will 
necessarily intersect a trap's region of influence. Now 
the least upper bound on the minimum value of |q — pm\ 
can be shown to be 1/{N + 1) Using the smallest 
allowed value of m, which is the slope of the line through 
the origin and tangent to the circle centered at (1,0), we 
estimate for an upper bound on Sm 



< 



1 



N+1 



(2) 



Thus a sufficient (and necessary) condition that any 
straight line trajectory starting at the origin must in- 
tersect a trap before exiting the lattice is Sm = a, or 



V > 



1 



(3) 



For the remainder of the discussion we shall assume that 
this condition is satisfied and that every trajectory in- 
volves a sequence of inter-trap jumps. 

Every trajectory starting at a lattice point on the left 
side of the lattice and exiting at the right side can be 
decomposed into a discrete sequence of steps from one 
lattice point to another, followed by exit from the trap 
array. For ease of modeling we treat the exit also as 
a step described by the lattice vector to the trap that 
would have captured the exiting particle had the trap 
array been extended to infinity in the x direction. 




tanG 

FIG. 2; Direction of deterministically locked-in transport, 
tan-f/), as a function of driving direction tanS. Results shown 
are for an array with trap density 77 = 0/6 = 0. 15. 



In the physical system |5| , colloidal particles are driven 
through an array of optical traps by a steady flow in the 
surrounding medium oriented at angle to the trap ar- 
ray's [10] axis. These particles are also subject to ran- 
dom thermal forces causing them to diffuse. We will first 
consider the case where there is no diffusion. In the ab- 
sence of traps, the particles then would travel in straight 
lines in the 9 direction. Within the trap array, however, 
a particle's displacement in each step is determined by 
the position of the center of the next trap it encounters 
along its straight line trajectory. Thus, even in the ab- 
sence of diffusion, the effective direction of motion can 
differ from the direction dictated by the flow. The effec- 
tive direction of motion will be described by an angle ^ 
given by tani/; = iy/ix where £ = {ix,£y) is the lattice 
vector describing the step. The functional dependence of 
ip on the flow direction, 9, is given by tan 7/; — q/p, where 
(p, q) are the smallest integers satisfying 



\q-pta.n9\ < 77 (1 + tan 6*^)5 



(4) 



with p E {I, - ■ ■ , N}. Figure 121 shows the result. 

Corresponding to a given value of 9 there exists a 
unique value of ip and hence a lattice vector, £, describ- 
ing that step. In general, a trajectory T consists of a se- 
quence of steps, T = {i^,£^, - ■ ■ , i^}, where ^ is a lattice 
vector describing the z-th step of an n-step trajectory. 
Here t = {P^,£l), with e {I,-- - ,iV} and (."y e Z, 
such that £^ and ty are relatively prime, ^ 
and X)r=i^ < These conditions ensure that we do 
not count multiples of lattice vectors as distinct and that 
each trajectory terminates at the right side of the trap 
array. We can now deflne a mean effective direction, ip, 
for the entire trajectory by 



tan7/>(r) 



V" pi 



(5) 



This value characterizes the particles' overall transport 
through the array, and would be reflected in measure- 
ments of bulk transport properties, such as a Hall voltage 
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in the case of a periodically modulated two-dimensional 
electron gas 0- Without diffusion or some other ran- 
domizing process, all the steps in every trajectory would 
be identical, with ijj = ip for any orientation 9. 

We now consider the case in which diffusion or scat- 
tering randomizes the particles' trajectories. The natu- 
ral assumption is that randomness would smear out the 
noise-free transport characteristic in Fig. [21 until, on av- 
erage, the trajectories simply follow the driving direction, 
ip = 9. In fact, far more interesting behavior results. 

We model trajectory fluctuations by assuming that a 
particle leaving a trap is equally likely to travel along any 
direction in a wedge of opening angle 2 S9 centered on the 
flow direction. In terms of experimental parameters for 
colloidal transport, 69 depends on the temperature and 
the particles' mobility. The path a particle takes through 
the array is no longer deterministic, and there are many 
possible trajectories for any given driving direction. 

We first compute the probability of a step being a lat- 
tice vector / = {ix,iy) with 4 G {1, • • • , A^} and £y € 
such that £x and iy are relatively prime. This proba- 
bility corresponds to the fraction of angles in the range 
[9 — 69, 9 + 69] whose effective direction ip corresponds to 
lattice vector £. The probability is given by 



2.0 F 



- ^9 



-50 



,tan-!/i(e') 



d9'. 



(6) 



where 6j,k is the Kronecker delta function. If the con- 
dition in Eq. is satisfied then there will only be a 
finite number of lattice vectors with non-zero probability 
as defined above, and the sum of the probabilities will 
be unity. Consequently there is only a finite number of 
possible trajectories consisting of sequences of these lat- 
tice vectors. Let there be M distinct possible trajectories 
given by the set {TJ. Here T, = {/i,^^, . . . ,/"■}, where 
rii is the number of steps in trajectory i. The probability 
for a particle to take a particular trajectory Ti is given 
by 



(7) 



The probability that the particle's trajectory will carry 
it in direction t/iqj is then given by 



M 



tan T/jQitan ipiTi) ' 



(8) 



This distribution is nonzero for only a few discrete direc- 
tions and thus is better characterized by its most prob- 
able value, V'*, with P{tp*) > P{ip) V ip, rather than its 
mean. 

Figure El shows the results for an = 9 lattice with 
ry = 0.26 when the dispersal angle 69 that models disor- 
der is increased from the deterministic limit, 69 = to 
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FIG. 3; Statistically locked-in transport through an A'^ = 9 
lattice with i] — 0.26 and 59 — 0.09 rad (circles), compared 
with deterministic transport (lines) and with experimental 
data from Fig. 3 of Ref. (squares). Data for statistically 
locked-in transport are offset by 0.5, and experimental data 
are offset by 1 for clarity. 



66 — 0.09. These values were selected to mimic the exper- 
imental conditions in Ref. j^. Data from that study are 
plotted as squares in Fig. O for comparison. The calcu- 
lated results demonstrate that disorder does not necessar- 
ily wipe out the structure of kinetically locked-in trans- 
port, but rather reconfigures the pattern of plateaus. In 
particular, not all statistically locked-in plateaus are cen- 
tered on commensurate directions. This was one of the 
principal surprises to emerge from experimental studies 
of colloid flowing through optical trap arrays and is 
a clear feature of the experimental data in Fig. |31 Now 
this restructuring can be explained as the emergence of 
high-order locked-in plateaus from a statistical sampling 
of low-order jumps. Indeed, as the bottom trace in Fig.|2| 
demonstrates, many of the commensurate directions cor- 
responding to statistically locked-in plateaus are not mi- 
croscopically accessible at 77 = 0.26. These plateaus are 
absent, therefore, from the underlying pattern of deter- 
ministically locked-in steps. 

By contrast to the higher-order states, the principal 
plateaus at q/p — and q/p — 1 are both microscopi- 
cally and deterministically locked-in. These correspond 
to transport along the [10] and [11] lattice directions, re- 
spectively, which cannot be decomposed into lower-index 
jumps. This distinction also is clear in the experimental 
data 0. Under some conditions, other directions such as 
[21] also can become deterministically locked-in, partic- 
ularly when both 77 and 69 are small. 

The calculated transport characteristics do not agree 
with the experimental data in all detail. Qualitative dif- 
ferences, such as the jump in the experimental data at 
q/p = 1/2 can be ascribed to the 10 x 10 structure of the 
experiment's optical trap array, which differs from the ex- 
tended geometry we considered here. Flowing colloidal 
particles also are not drawn directly to optical traps' cen- 
ters as they jump through a trap array, but rather fol- 
low more subtle and complicated trajectories. Clearly, 
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FIG. 4: Dependence of the most probable transport direction, 
tan ■(/)*, on driving direction ta.nd in the presence of random 
forcing. Circles show results for with = 0.15, 59 — 7r/24 
and A'^ = 10. The solid curve shows equivalent results for 
ri — 0.4 and 59 = it/32, and illustrates the extent to which 
qualitative features of the transport can change with the size 
of the flowing objects. 

though these measurements more closely resemble our 
model's predictions than they resemble the simple deter- 
ministic spectrum for the same conditions. 

The distribution of statistically locked-in states 
changes dramatically with the size of the array. Acces- 
sible states proliferate as N increases, first smearing out 
the spectrum of statistically locked-in states, and then 
erasing even the deterministic plateaus. In the large-A^ 
limit and for 66 > 0, the array has no overall effect on 
transport: = %p = 6. The hierarchically structured 
transport characteristics discovered experimentally and 
explained here are features of mesoscopic systems. 

Being able to tune transport properties by adjusting 
the number, distribution, and size of discrete potential 
wells suggests a practical method for continuously sep- 
arating heterogeneous fluid-borne samples into multiple 
fractions simultaneously. Consider, for example, objects 
driven through arrays of discrete optical tweezers 0, 
each of which is formed by a strongly focused beam of 
light. Different objects passing through this optical in- 
tensity field experience different potential energy land- 
scapes depending on their size, shape, and composition. 
This effect already has been exploited to sort mixtures 
into two fractions with a single line of optical traps [loj . 
In terms of our model, the disk-like wells would appear 



to be larger for larger or more strongly trapped objects. 
Larger objects also would diffuse less vigorously, and so 
would sample a simpler spectrum of statistically locked- 
in states. Consequently, the same array of potential wells 
can have a substantially different effect on different ob- 
jects' transport. 

Figure^ldemonstrates that smaller particles can be de- 
flected into high-order statistically locked in states under 
the same conditions that larger objects are determinis- 
tically locked in to the [10] direction. The resulting an- 
gular separation is large enough that the two resulting 
fractions could be collected, for example in microfluidic 
channels. Indeed, this figure makes clear that several 
fractions might be continuously selected out of a mixed 
sample by a single static array of optical traps. 
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